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ABSTRACT 



Context. The presence of water in the wind of the extreme carbon star IRC+10216 has been confirmed by the Herschel telescope. 
The regions where the high-J H2O lines have been detected are close to the star at radii r < 15 Ri,. 

Aims. We investigate the formation of water and related molecules in the periodically-shocked inner layers of IRC+10216 where dust 
also forms and accelerates the wind. 

Methods. We describe the molecular formation by a chemical kinetic network involving carbon-and oxygen-based molecules. We 
then apply this network to the physical conditions pertaining to the dust-formation zone which experiences the passage of pulsation- 
driven shocks between 1 and 5 R^. We solve for a system of stiff, coupled, ordinary, and differential equations. 
Results. Non-equilibrium chemistry prevails in the dust-formation zone. H2O forms quickly above the photosphere from the synthesis 
of hydroxyl OH induced by the thermal fragmentation of CO in the hot post-shock gas. The derived abundance with respect to H2 at 
5 i?* is 1.4x10"^, which excellently agrees the values derived from Herschel observations. The non-equilibrium formation process of 
water will be active whatever the stellar C/O ratio, and H2O should then be present in the wind acceleration zone of all stars on the 
Asymptotic Giant Branch. 
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1. Introduction 

' The extreme carbon star IRC+10216 is one of the best studied 
, evolved low-mass stars on the Asymptotic Giant Branch (AGB) 
' owing to its proximity (d ~ 180 pc) and the rich chemistry of 
its wind. Indeed, more than 60 molecular species have been de- 
tected at millimetre (mm), submillimetre (submm) and infrared 
(IR) wavelengths, probing the entire gas conditions pertaining to 
the outflow (Ziurvs l2006l l. Because the star has already experi- 
enced third dredge-up, it is 'carbon-rich', i.e., its photosphere 
is characterised by a C/O ratio greater than 1. Many carbon- 
bearing species have been identified in the inner envelope ex- 
tending from 1 to ~ 10 Ri,, including CO, HCN, C2H2, CS, 
and SiC2 at mid-IR and submm wavelengths (Keady & Rigdway 
!ri993 Schoier et al. 2007, Fonfria et al. 2008, Patel et al. |2009[ 
' Cernicharo et al.^OlO). Oxygen-bearing species other than CO 
I have also long been detected in IRC+10216. Submm observa- 
. tions of silicon monoxide, SiO, confirmed a formation locus 
close to the star (Schoier et al. 2006, Decin et al. 120101 ). The wa- 
ter molecule, H2O, was first detected with the submm satellite 
SWAS by Melnick et al. (2001), and its presence was confirmed 
by submm observations with the ODIN satellite (Hasegawa et al. 
I2OO6I 1 and supported by the detection of hydroxyl, OH, by Ford 
et al. (2004). Recently, water was detected by the SPIRE, PACS, 
and HIFI spectrometers onboard the Herschel submm telescope 
in the carbon star VY Cyg (Neufeld et al.l2010K the S star^ Cyg 
(characterised by a C/O ratio ~ 1) (Justtanont et al. 120101) . and 
IRC+10216 (Decin et al. 2010, hereafter DABIO). These vai'i- 
ous detections of high-excitation rotational lines probe relatively 
high gas temperatures, thus implying the presence of H2O fairly 
close to the star. Derived abundances with respect to H2 span 



values from 10"** to 10"^ for IRC+10216, ~ 10"^ for VY Cyg, 
and 10"^ for;^' Cyg. The derived abundances seem to logically 
decrease as the carbon content of the star increases, because CO 
traps most of the oxygen available for water formation. 

Several formation mechanisms were advocated to explain 
the presence of H2O in IRC+10216. Melnick et al. (12001 ) pro- 
posed that icy comet bodies orbiting the carbon star could be 
vaporised in the stellar outflow, providing a source of oxygen 
to the gas for water formation in the intermediate envelope. 
Willacy (2004 ) suggested that water could form on the surface of 
iron dust grains in the intermediate envelope by Fischer- Tropsch 
catalysis. Recently, DABIO and Agundez et al. ("2010 hereafter 
ACGIO) proposed that partial penetration of the interstellar ul- 
traviolet (UV) radiation field occurs deep in the outflow owing to 
the clumpy nature of the wind. '^CO and SiO can thus photodis- 
sociate, providing atomic oxygen to the gas, which then leads 
to the formation of H2O. All the proposed explanations for the 
presence of water have their drawbacks. The first mechanism is 
somehow extreme because it implies that orbiting icy cometary 
bodies should be a characteristic of all carbon and S stars. The 
second proposition requires that iron grains form in the wind of 
carbon stars, an assumption that still needs confirmation by the- 
oretical models or observations. Finally, according to the third 
model, all stellar winds are sufficiently clumpy to allow for some 
penetration of UV photons as deep as 2 7?*, that is, in the dust 
formation and wind acceleration region. If so, partial photodis- 
sociation of molecular dust precursors, among which the radical 
propargyl (C3H3) and its precursors CH and CH2, should occur, 
which would hamper the dust-formation process to a certain de- 
gree. 
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Table 1. Gas parameters in the inner wind of IRC+10216. The shocks form at 1.2 7?*. For each radius, the temperature and number 
density are given for the pre-shock gas, the gas after the shock front (in the collision-induced H2 dissociation zone), and the gas at 
the beginning and the end of the adiabatic cooling zone (= ballistic trajectory). 



r 


Shock velocity 


Preshock gas 


Shock Front 


Adiabatic 


expansion - start 


Adiabatic 


expansion - end 




(kms-') 


T 




T 


^gas 


T 




T 




1.2 


20.0 


2062 


3.63(13) 


19725 


1.98(14) 


4409 


5.97(14) 


1480 


3.63(13) 


1.5 


17.9 


1803 


8.24(12) 


15922 


4.40(13) 


3870 


1.29(14) 


1290 


8.24(12) 


2 


15.5 


1517 


1.44(12) 


12081 


7.59(12) 


3200 


2.14(13) 


1080 


1.44(13) 


2.5 


13.9 


1327 


4.24(11) 


9779 


2.21(12) 


2750 


6.08(12) 


951 


4.24(11) 


3 


12.6 


1190 


1.69(11) 


8245 


8.73(11) 


2430 


2.35(12) 


848 


1.69(11) 


4 


11.0 


1001 


4.51(10) 


6284 


2.29(11) 


1790 


4.48(11) 


711 


4.51(10) 


5 


9.8 


876 


1.79(10) 


5096 


8.94(10) 


1550 


1.71(11) 


621 


1.79(10) 



Notes. Temperatures T are in Kelvin and gas number densities n^a, are in cm ^ . 



The formation of the oxygen-bearing species SiO in the inner 
wind of IRC -1-12016 was investigated theoretically by Willacy & 
Cherchneff (119981 hereafter WC98), who showed that coUisional 
dissociation of CO occurred in the gas layers that experienced 
the passage of periodic shocks induced by stellar pulsation. The 
released atomic oxygen then formed a population of OH radi- 
cals that triggered the synthesis of SiO via their reaction with 
atomic Si. Shock-induced chemistry could also explain the for- 
mation of CO2 in the O-rich Mira star, IK Tau (Duari et al. ll999l) . 
In a later study, Cherchneff (2006) modelled the non-equilibrium 
chemistry of the inner wind of AGB stars as a function of C/O ra- 
tios. Of importance was the finding that a few molecules, namely 
CO, SiO, HCN, and CS, efficiently formed in large amounts in 
the dust-formation zone whatever the C/O ratio of the star. The 
author concluded that this group of molecules was ejected as 
parent species in the intermediate and outer envelopes in O-rich 
Miras, S stars, and carbon stars. This hypothesis was additionally 
supported by the observations of high-excitation rotational lines 
of CO, HCN, SiO, and CS in O-rich, C-rich, and S AGB stars 
(Decin et al. 2008). These results strongly support the assump- 
tion that the non-equilibrium chemistry induced by the passage 
of periodic shocks is responsible for the formation of C-bearing 
species in the inner envelope of O-rich Miras, and of O-bearing 
species in the inner wind of carbon stars. 

Here, we revisit the non-equilibirum chemistry of the inner 
wind of IRC-i-10216. The updated chemistry includes new pro- 
cesses such as the thermal fragmentation that is active at the high 
postshock gas temperatures and the formation and destruction 
of O-bearing and C-bearing species. The latter include acetylene 
C2H2, hydrocarbons, carbon chains, and the benzene and phenyl 
aromatic rings, CgHe and CeHs, respectively. The chemistry 
of metal hydrides, chlorides, and sulphides, and phosphorous- 
bearing compounds is also considered. The complete results 
of this study will be presented in a forthcoming publication, 
and we report here on our results for water. Section 2 presents 
the physical and chemical model considered for the inner wind 
of IRC+10216, while our results for H2O and other important 
species are summarised in Section 3, and a discussion is pre- 
sented in Section 4. 

2. The physical and chemical model 

The dust-formation zone of IRC+ 12016 extends above the pho- 
tosphere from 1 to ~ 5 Ri,, according to Keady & Ridgway 
(II993I 1. For the sake of clarity, we refer to this region as the 
'inner wind' in the rest of the paper. These layers experience 
the periodic passage of shocks induced by the stellar pulsation. 
Following Bowen ( I1988I I. these shocks compress the gas lay- 



Table 2. Chemical species that are linked to the water chemistry 
and are included in the inner wind model of IRC+10216. 



Atoms 


H 





C 


Si 


S 


N 


Diatomic 


H, 


CO 


SiO 


OH 


O2 


SO 




NO 


C2 


CS 


CN 


CH 


SiS 




S2 


N2 


NH 








Tri-atomic 


H2O 


HCN 


CH2 


C2H 


HCO 


C3 




CO2 


NH2 










4-atoms 


CH3 


NH3 


C2H2 








Hydrocarbons 


C2H3 


C3H2 


C3H3 


C4H2 


C4H3 


C4H4 


Aromatics 




CeHe 











ers, resulting in high temperatures and densities, which further 
relax though the collisional dissociation of H2 (Fox & Wood 
119851) and adiabatic expansion (Bertschinger & Chevalier ll985l l. 
Experiencing the stellar gravitational field, these gas layers fall 
back to almost their initial position. We model the inner wind 
following the formalism used by WC98 and Cherchneff (2006jl, 
who studied the chemistry of the immediate postshock gas cool- 
ing by H2 dissociation and subsequent adiabatic expansion. We 
assume that the shocks form at a radius rj = 1.2 above 
the photosphere with an initial strength of 20 km s ' (Ridgway 
& Keady I1981II. Th e C/O ratio is taken to be equal to 1.4 
(Winters et al. 119941). The initial pre-shock gas density at the 
shock-formation radius was rescaled down by a factor of ten 
(ngasifs) - 3.63 X lO'-' cm""*) compared to the initial value used 
by WC98, which was considered to be too high (Agilndez et al. 
2006). We will see below that this change has no impact on the 
formed chemical species and their respective abundances. The 
gas is followed over one pulsation period P from 1 to 5 and 
the preshock and postshock gas parameters as a function of ra- 
dius are listed in Table [1] 

We updated the chemistry of the inner wind and included all 
chemical processes relevant to the high temperatures and num- 
ber densities characteristic of the shocked gas layers. All chem- 
ical pathways leading to the formation of linear molecules, car- 
bon chains, and aromatic rings, include neutral-neutral processes 
like termolecular, bimolecular reactions and radiative associa- 
tion reactions, whereas destruction is described by thermal frag- 
mentation and neutral-neutral processes (i.e., oxidation reactions 
of hydrocarbons and all reverse processes of the formation re- 
actions). No ions are considered in this chemistry because the 
UV stellar radiation field of IRC+10216 is too low to foster ef- 
ficient photodissociation and ionisation processes. The periodic 
shocks also have insufficiently high velocities to be radiative, 
and as previously mentioned, the immediate postshock gas cools 
via H2 collisional dissociation (Fox & Wood J985) . The rates 
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Fig. 1. H2O abundances with respect to total gas as a function of 
pulsation phase and radius in the inner wind (shocks form at 
61 = and 61= 1). 
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Fig. 2. Abundances with respect to H2 of key molecules includ- 
ing H2O taken at 9 - 1, as a function of radius. Values at r = 1 
R-i, are derived from thermodynamical equilibrium calculations. 



for the chemical processes are taken from the National Institute 
of Standards and Technology database (NIST chemical kinetics 
database) and also from the literature of combustion, atmopsh- 
eric, and material sciences. Species linked to the shock chem- 
istry of H2O in the inner wind are listed in Table |2] and more 
detail on the chemistry is given in the online Appendix A. 

3. Results 

The abundances of H2O with respect to total gas are presented 
in Figure [T] as a function of radius in the inner wind and pul- 
sation phase - tjP where t is the time and P the pulsation 
period. The results follow the trends highlighted by our previous 
studies (WC98, Cherchneff 2006). A non-equilibirum chemistry 
takes place in the post-shock gas at radius 1.2 resulting in 
the coUisional breaking of CO at the high post-shock tempera- 
tures. H2 is also destroyed in the immediate post-shock gas and 
reforms in the adiabatic expansion at phase > 0.3. The released 
atomic oxygen reacts with H2 to form hydroxyl according to the 
reaction 

+ H2—^0H + H, (1) 

whose rate has an energy barrier of ~ 3200 K. Therefore, this 
reaction is efficient at the high post-shock gas temperatures en- 
countered close to the star (see Table [1]). Hydroxyl additionally 
reacts with H2 to form H2O according to 

0H + H2-^ H2O + H. (2) 

Reaction |2] has a lower energy barrier and a faster rate than 
Reaction [Hand allows for H2O formation at lower temperatures. 
Close to the star, the hot post-shock gas temperatures preclude 
the reformation of H2 at early phases and the efficient formation 
of H2O according to the above processes, leading to very small 
water abundances. H2O forms at > 0.4 once H2 has reformed. 
At radii exceeding 4 the shock velocity is low and such 
that almost no destruction of H2 occurs at early phases, which 
triggers the early processing of H2O. Once formed, H2O abun- 
dances remain constant as the low gas temperatures freeze out 
the neutral-neutral, high-T chemistry. In Figure|2]the abundances 
with respect to H2 of specific species are shown as a function of 
radius in the inner wind. These molecules (CO, HCN, SiO, and 
CS) belong to the group of species predicted to always form with 



large abundances in the inner wind of stars during their evolution 
from the Mira to the carbon star stage (Cherchneff 2006). The 
abundance of H2O shows similar trends as that of SiO because 
those two molecules are competitors in the consumption of OH 
and their chemistry is coupled (see Appendix A). 

A striking result is the excellent agreement between our pre- 
dicted water abundances at radii > 3 Ri, and those derived by 
Melnick et al. (2001) from their SWAS observations and by 
DABIO from their Herscliel/SPIRE/PACS data. At 5 we de- 
rive an abundance with respect to H2 of 1 .4 x 10"^ where the val- 
ues from the SWAS and the Herschel/HIFI data are ~ 1x10"^. 
In previous studies (WC98, Cherchneff I2006I I. H2O abundances 
showed similar trends as the present ones, i.e. a very efficient 
formation at the shock-formation radius and a decrease farther 
away from the star. However, the decline at larger radii was 
much sharper than here. This is essentially because of the greater 
completeness of the chemistry considered in the present model. 
A more accurate chemistry was considered for silicon and sul- 
phur, leading to lower abundances of SiO than previously found. 
These SiO abundances with respect to H2 level off at a value of 
~ 6 X 10"*^ at r = 5 and perfectly agree with the value range of 
2 X 10"^ - 3 X 10"^ derived from the new Herschel/PACs data on 
Si-bearing species (Decin et al. 2010 ). The lower SiO content re- 
sults in more available OH radicals to form water, as described in 
Appendix A. The formation of the aromatic molecule benzene, 
Cdlid, and its radical phenyl, CeHs, as well as their hydrocar- 
bon precursors is specifically taken into account along with their 
destruction by oxidation reactions. These destruction processes 
involve OH, the prevalent oxidation agent in the wind, and result 
in the formation of water. Finally, in carbon stars, a great fraction 
of the carbon not locked in CO goes into hydrocarbons, starting 
with acetylene, C2H2, whose calculated abundances is slightly 
lower than 1 x 10 "* in the inner wind. The formation and de- 
struction of C2H2 and larger hydrocarbons impacts on the abun- 
dance of molecular hydrogen H2 and indirectly on the formation 
of water, as explained in Appendix A. These combined effects 
coupled to the freeze-out of H2O chemistry at r > 3 7?* lead to 
higher H2O abundance values than in previous studies at these 
radii. 
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4. Discussion 

The models of ACGIO and DAB 10 allow for the penetration of 
the interstellar UV radiation field in the very deep layers of the 
inner wind owing to clumping in the outflow. The models use 
that of Agundez & Cernicharo ( I2006I I. i.e., a Eulerian descrip- 
tion of a steady outflow penetrated by some UV radiation (in the 
form of a minor UV-illuminated wind component), and do not 
take into account the periodic shocks pervading the inner wind 
from 1 to ~ 5 Ri,. The region exposed to UV radiation is 10 
% of the stellar envelope mass, and the models are highly de- 
pendent on the clumping factor and stellar mass-loss. The free 
oxygen necessary for water synthesis is provided by the partial 
photodissociation of '^'CO and SiO. This mechanism also trig- 
gers the formation of other species of interest, like ammonia. 
Indeed, NH3 is predicted to form at radius r > 3 ^* with high 
abundances peaking at 7 x 10"^ with respect to H2 at r ~ 4 
Keady & Ridgway (119931 observed mid-IR vibrational transi- 
tions of ammonia in IRC+10216 and deduced that a NH3 abun- 
dance distribution peaking at 10-20 Ri, could better reproduce 
their data. Later, Monnier et al. (120001) carried out interferomet- 
ric observations of mid-IR molecular absorption bands of NH3 
in IRC+10216 with very high spectral resolution. They found 
that NH3 was located in a region of decaying gas turbulence at 
radii well beyond the inner wind (r > 20 Ri,). These large radii 
were further used to model NH3 submm lines observed with the 
satellite ODIN (Hagesawa et al. 120061) . The above-mentioned 
observations suggest that the formation locus of NH3 is located 
well beyond the dust-formation zone (r > 5 Ri,), while ammonia 
is formed as early as 3 /?*with a high abundance extending to 
~100 according to DAB 10 and ACGIO. Our present model 
does not form NH3 in the inner wind (ji:(NH3) ~ 4 x 10 with 
respect to H2 at r = 5 R*), supporting larger formation radii for 
this species in accordance with mid-IR and submm observations. 

Decin et al. ( 120101 DAB 10) report on newly observed high- 
excitation lines of cyanoacetylene HC3N with the IRAM tele- 
scope. These high-J lines present a flat-topped profile, and 
DAB 10 and ACGIO construe this shape as evidence for a for- 
mation locus deeper in the envelope. They claim that the dis- 
sociation occurring in their minor UV-illuminated wind compo- 
nent explains the formation of HC3N in the intermediate enve- 
lope. An enhancement (up to ~ 3 x 10"^) in the HC3N abun- 
dance peaking at r = 120 is required to reproduce the high-J 
IRAM lines. The radial distribution of HC3N in IRC+10216 was 
mapped by Audinos et al. (119941 ) with the IRAM telescope, and 
the line modelling for a steady homogeneous outflow already re- 
quired an enhanced HC3N abundance at these radii. Cherchneff' 
& Glassgold (1993) modelled the cyanopolyyne chemistry in 
IRC+10216 and found a shoulder in the HC3N distribution at 
radii < 200 Ri, resulting from the synthesis of HC3N by neutral- 
neutral channels. The predicted HC3N abundances in the shoul- 
der were a factor ~ 7 lower compared to values derived by 
Audinos et al. We notice that the HC3N abundance distribu- 
tion of DAB 10 for their major UV-shielded component does 
not include this predicted shoulder, which points to a different 
chemistry used by DAB 10 and ACGIO. High angular resolution 
observations by Trung & Lim (2008) of the J = 5 - 4 HC3N 
line with the VLA clearly indicate a very clumpy shell distri- 
bution as far as 20" from the star. Therefore, radiative transfer 
models taking into account the outflow dumpiness are required 
to accurately assess the HC3N abundance distribution across 
the envelope. An updated chemical model is also necessary to 
clearly identify the prevalent formation pathways to HC3N, in- 
cluding neutral-neutral and UV dissociation-induced processes. 



A shock-induced chemistry may also be a viable source of HC3N 
at intermediate radii and must be considered in subsequent mod- 
els. 

The present results support the hypothesis that H2O forms 
very close to the star with high abundances with respect to H2 
between 1 x 10"*' and 1 x 10""* and that the abundance grad- 
ually chemically freezes out to a value of 1.4 x 10"^ at r > 5 
Ri,. The contribution to the H2O line intensities from the high- 
abundance region comprised between 1.2 Ri, and 2.5 Ri, is not 
observable with Herschel because of its small extent and beam 
filling factor. Furthermore, the H2O formation chemistry appears 
to be coupled to that of SiO and highlights the importance of 
the hydroxyl radical OH in the inner wind. A major advan- 
tage of the shock-induced chemistry hypothesis is that any star 
on the AGB pulsates and experiences the passage of shocks in 
its dust-formation zone. Therefore, it provides a non-restrictive, 
genuine mechanism for forming water and other molecules very 
close to the star for a variety of objects in different evolution- 
ary stages, potentially explaining the detection of H2O in O- 
rich, S and carbon stars with the Herschel telescope. The shock- 
induced chemistry is also successful in explaining the presence 
of O-bearing molecules in carbon stars and C-bearing species in 
O-rich Miras. Foremost, this hypothesis does not exclude other 
formation mechanisms for water at larger radii (comets, chem- 
istry on grain surfaces, partial photodissociation in a clumpy 
outflowlJ Because of its presence in AGB stars, H2O can be 
added to the fist of species (CO, SiO, HCN, CS) proposed by 
Cherchneff (120061) to efficiently form in the dust-formation zone 
from shock chemistry and be ejected as 'parent' species in the 
intermediate and outer envelopes. This hypothesis awaits testing 
by observations of the deep layers of AGB envelopes at mid-IR 
and submm wavelengths. 
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Table A.l. Chemical reaction types included in the chemical 
model of IRC+10216 inner wind. 



Unimolecular AB — > A + B Decomposition 

Bimolecular A + B — > C + D Neutral-neutral 

A -I- B — > AB + hv Radiative association 

AB + M — > A -I- B + M Thermal fragmentation 

Termolecular A + B + M — > AB + M Three-body association 



Appendix A: Chemical model and new key 
processes 

The present chemical model is based on that of Willacy & 
Cherchneff (119981 hereafter WC98) and includes all processes 
relevant to the hot and dense gas pertaining to the inner wind be- 
tween 1 and 5 However, many processes and their rates have 
been updated and/or added in view of the progresses made in 
chemical kinetics, and combustion and aerosol chemistry since 
1998. The various types of chemical pathways considered in the 
model are summarised in Table lA. fl Thermal fragmentation, i.e., 
destruction of molecules by collision with the ambient gas, uni- 
molecular decomposition of aromatics and hydrocarbons, and 
radiative association reactions were added to the chemistry. The 
chemical network includes 59 species (some of which are listed 
in Table|2]i and 370 chemical reactions. The rates for these chem- 
ical processes are taken from the National Institute of Standards 
and Technology database (NIST chemical kinetics database), 
and from the literature of combustion, atmopsheric, and mate- 
rial sciences. 

In terms of formalism, three major changes were imple- 
mented with respect to WC98 and Cherchneff (I2006I I. Firstly, 
the treatment of the reverse reaction of a specific process was 
changed. Several new rates were measured in combustion and 
aerosol chemistry and are now available. Therefore, instead of 
calculating the rate of the reverse process from detailed balance 
and the equilibrium constant (see Equation 4 in WC98), we di- 
rectly enter the available measured or calculated rate values in 
the chemical network. When the information is not available, we 
make 'educated' guesses depending on the type of the reaction. 

Secondly, the chemistry involving atomic silicon. Si, and Si- 
bearing species has been restricted to reactions for which rates 
were measured or calculated. The reactions and rates derived 
from the isovalence of Si with carbon as stated in WC98 were 
abandoned. The prevalent formation reactions for SiO are 

Si + CO SiO + C, (A.l) 

and 

Si + OH SiO + H, (A.2) 

Reaction lA.ll is effective at forming SiO at temperatures > 
3000 K, whereas Reaction IA.2I contributes to most of the SiO 
synthesis at temperatures < 3000 K. Destruction of SiO is mainly 
triggered by the opposite process of Reaction IA.2l at all temper- 
atures. Reaction IA.2l is a competitive channel to the formation 
of H2O via Reaction|2] because it depletes OH radicals. But be- 
cause of the high H2 content of the wind, the net formation rate 
of H2O according to Reaction |2] is much higher than that for 
SiO. Furthermore, a large quantity of atomic Si is tied up in sili- 
con sulphide, SiS, leaving a reduced pool of available Si to form 
SiO. These combined effects result in water abundances higher 
than those of SiO, as seen in Figure |2] 

Thirdly, the chemistry now included the formation of a larger 
set of molecules, some of them listed in Table |2] As explained 



in Section 3, the formation of single aromatic ring compounds 
(benzene and phenyl) was considered. Previous studies of the 
formation of polycyclic aromatic hydrocarbons in AGB winds 
never took into account the presence of O-bearing species in 
the shocked regions (CherchneflF 120 1 Oi l . The coupling of these 
C- and O-rich chemistries is important, because the unimolecu- 
lar decomposition of hydrocarbons and aromatics replenish the 
gas in acetylene. Since C2H2 has large abundances in the inner 
wind and its formation and destruction processes are coupled to 
the H2 chemistry, it indirectly impacts on the formation of OH, 
H2O, and SiO through Reactions [1] |2] and IA.2I Furthermore, 
the potential oxidation of hydrocarbon precursors by O-bearing 
species, namely OH, is included and gives back H2O to the gas 
phase, albeit its impact is minor in replenishing water 



